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Abstract: Recently, Freyhult, Rej and Staudacher (FRS) proposed an integral equation 
determining the leading logarithmic term of the anomalous dimension of 5l{2) twist- 
operators in ^ = 4 SYM for large Lorentz spin M and twist L at fixed j = L/\og M. 
We discuss the large j limit of the FRS equation. This limit can be matched with the fast 
long string limit of AdS^ x superstring perturbation theory at all couplings. In par- 
ticular, a certain part of the classical and one-loop string result is known to be protected 
and can be computed in the weakly coupled large-^' limit of the FRS equation. We present 
various analytical and numerical results supporting agreement at one and two loops in 
the gauge theory. 



Contents 

|l[ Introduction Q 

^ The large-j limit and the fast long string limit ^ 

^ The FRS equation in brief § 

^ Setup I 

^ The all-loop FRS equation 



The two-loop hole density equation in u-space g 

H NNLO large-; expansion of the one-loop FRS equation g 

^3] Leading order ^ 

5!2| Next-to-leading order ^ 



5.3 Next-to-next-to-leading order 



|6|, Large-; expansion of the two loops FRS equation: LO and remarks 15 

1^ Numerical study of the hole density equation 16 

^ One loop y 



7.2| Two loops |17 
|8[ Conclusions 17 
^ Technical results about the combinations ip{a + ix) + ■tp{a — ix) 18 



1. Introduction 

The anomalous dimensions of Wilson twist operators ^ are relevant perturbative quan- 
tities which appear in various phenomenological problems in the study of QCD strong 
interactions. A typical example is the operator product expansion analysis of deep in- 
elastic scattering [Q]. In that context, the close relation between parton splitting functions 
and anomalous dimensions suggests various physical insights valid in special kinemati- 
cal limits. In particular, the behavior of anomalous dimensions for large Lorentz spin M 
at fixed twist L probes the quasi-elastic limit where the Bjorken variable is close to unity 
XBj 1. In this regime the most singular part of the splitting functions is due to soft gluon 
emission and is universal. For the leading twist 2 operators, these remarks translate into 
the following well-known prediction for the anomalous dimension 7 



7 = 2rcusp(5) log M + 0(M°), 



(1.1) 



where = and X = Nc 5ym Hooft planar coupling A. The non trivial function 

Tcuspl^) is the so-called cusp anomalous dimension [^]. 

The logarithmic scaling in Eq. is quite general and applies in particular to the 
superconformal finite JV = 4 SYM theory where integrability [^] and AdS/ CFT duality 
^ ^ can be exploited to gain (much) additional information. This approach is clearly 
interesting in itself due to the theoretical relevance of the J\A = 4 SYM theory. Besides, 
one can also argue that large XBj physics can be related to the QCD one, being mostly 
related to the shared gauge sector. A recent example of this strategy is the analysis of a 
generalized Gribov-Lipatov reciprocity ^ |o|] for various twist-2 and twist-3 Wilson 
operators [0, ||, |l|]. 

The current knowledge of rcusp(5) in A/" = 4 SYM is quite complete. It can be ex- 
tracted from the anomalous dimensions of 5l(2) operators. The weak-coupling perturba- 
tive series can be computed at all-orders by a rather simple expansion of the so-called BES 
equation [|T6|]. The result is in agreement with the most advanced available field theoret- 
ical computations |17]. The problem of computing the strong coupling expansion of the 
BES equation is more difficult and after intense activity [18] has been impressively solved 
in the remarkable paper [[T9[]. Again, there is full agreement with the two-loop analysis of 



the dual superstring theory on AdSs x [gO], gl], 

The BES equation is derived by considering operators with an arbitrary finite twist 
L and taking the large spin limit M — ^ co. If the twist increases with the spin M then 
one expects a richer landscape of scaling behaviors. A simple one-loop illustration of this 
general statement can be found in [24]. It is shown that when M ^ L, one must still 
distinguish between two quite different regimes characterized by extreme values of the 
gauge theory parameter ^ defined as 

In particular, the minimal anomalous dimension has the following leading contributions 

r log M, 

The first case is covered by the BES equation. The second case with the characteristic 
double logarithm enhancement is beyond its reach. The appearance of these two regimes 
is in quite similarity with the semiclassical string calculation of [ pi] ] as we shell discuss in 
a moment. 

In [^, Freyhult, Rej and Staudacher (FRS) proposed to analyze the logarithmic be- 
havior of anomalous dimensions in the following limit 

L.M ^ 00, j = — = fixed. (1.4) 
log M 

In this limit FRS prove that a logarithmic scaling is observed once more. The prefactor 
now depends on both g and j 

7(5,J)=/(5,J) log M + 0(M°), (1.5) 
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where f{g,j) is a generalization of the cusp anomalous dimension 



/(5,0) =/(5) = 2rcusp(5)- (1-6) 
An integral equation analogous to the BES equation, but valid for all g and j, has been 



derived in [p5[]. Of course, a great deal of interesting results can be obtained by applying 
to the FRS equation the methods which have been already sharpened in the case j = 0. 
In particular, this means that the FRS equation can be considered in the following two 
opposite limits. 

1. The fully weak limit. This is simply g,j 0. There seems to be no ambiguity in this 
double limit and it is convenient to first expand f{g,j) around j = 

f[g,i) = Y.fr.{g)r, (1-7) 

and then expand each fn{g) around 5 = 

/n(5) = E-^n,fc5''- (1-8) 

fc>0 

The coefficients Tn,k have been computed in p5| ] where their explicit expression can be 
found as well as a discussion of various features, like for instance transcendentality uni- 
formity. 

2. The Alday-Maldacena limit. The Alday-Maldacena (AM) limit is a strong coupling limit 
defined by the general condition g oo with j <^ g [^. The scaling function f{g,j) is 
described in this limit by the thermodynamical Bethe Ansatz equations of the non linear 
0(6) (T-model [^. As explained in the beautiful analysis of [^] it is necessary to con- 
sider separately the two situations where j <^ m or j ^ m where m is the dynamically 
generated mass gap [ ^ ] 

o3/4 1/4 

In particular, the case j m <^ g predicts the large g behavior of the functions fn{g) and 
can be summarized by the expansion 



fig,j) = -J + m- 



2 



m 24 Vm 



(1.10) 



which has been indeed recovered in the FRS equation in pq]. For additional numerical 



and analytical confirmations of the expansion Eq. ( I.IOD see |3l]]. Additional terms in 



the above series which represent the a-model energy density can be found in [^]. The 
other limit m ^ ^ 5 is also very interesting and is discussed in details in [^]. 

The above two limits are similar to those already considered for the cusp anomalous 
dimension since the parameter j is used as a perturbative book-keeping device. This 
suggests to consider another new limit. 
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3. The large- j limit. A quite different and very interesting limit is obtained taking first the 
weak coupling perturbative expansion of f{g,j) around 5 = 



/(5,j) = E/'"^(»5'" (1-11) 

n>0 

The functions f^'^\j) can be expanded around j = recovering the fully weak regime. 
On the other hand, one can look for the large j behavior of /'"^ (j). This large-j limit turns 
out to be non trivial. Looking back at the analysis of p^], we see that for large M and 
fixed j we simply have 

^ = -. (1.12) 
J 

The result Eq. ( |1.3D can be nicely rewritten in a uniform way as 

«1) = 8, /(!)(; »i) = ^. (1.13) 

J 

These simple relations immediately suggest that the large-j limit of the FRS equation is 
closely connected to the string theory calculations described first in [ |2l] ] and later ex- 
panded in [33]. In particular, the string perturbative calculations admit a BMN-like ex- 
pansion which is captured by the FRS equation in the large-j limit. The comparison can 
be done at arbitrary coupling, thus going beyond Eq. ( [1.11 ). This is important if one is in- 



terested in detecting universal dressing effects. In this paper, we exploit the fact that the 
BMN-like expansion contains some terms which are protected and can be computed in 
the weakly coupled gauge theory. This means that they can be matched by studying the 
large-j expansion of the one and two-loop expansion of the FRS equation, a remarkable 
simplification. Thus, we analyze the large-j limit by analytical and numerical methods 
and provide strong support for some of the predictions following from the computations 

The plan of the work is the following. In Sec. we briefly recall a few basic facts 
about the so-called fast long string limit of the folded string solution. In Sec. (^) we give 
a self-contained summary of the FRS equation. In Sec. (^) we present the explicit two- 
loop hole density equation in Bethe roots space. In Sec. (^) and Sec. (§) we describe the 
analytical re-derivation of the one-loop large-ji limit at next-to-next-to-leading order and 
a few considerations about the similar analysis at two loops. Finally, in Sec. @ we show 
our numerical results supporting at two loops the predictions of [|33|]. 
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2. The large-j limit and the fast long string limit 



In [21, S. Frolov and A. A. Tseytlin compute the semiclassical expansion around the 
rotating folded string configuration extending the analysis of [20, 34 1 and including the 
string center of mass motion along a big circle of S^. Their solution depends on the 
Lorentz spin M and 50(6) spin L to be identified with the quantum numbers of the s[(2) 
twist operators. The large A expansion of the energy takes the usual form 



Bo 

For an alternative derivation of the one-loop contribution see also [^]. The expansion 
Eq. ( [2.1[ ) can be considered in the long string limit which is 

long string both 1 < ^ and ^ < (2.2) 

vA vA vA 

This limit can be further refined in the two sub-cases defined by the additional conditions 

L M 

slow long string C log (2.3) 
v A vA 

fast long string log^^C-7=. ^^'^^ 
V A V A 

An interpolating regime between these two cases is obtained by fixing the parameter 

x = — log—. (2.5) 

TT L L 

The slow long string limit is reproduced for a; ^ 1, the fast string limit for a; <^ 1. The 
expansion of the energy in this second limit reads [|^] (see [36| for a recent analysis of E2 
mainly at a; C 1) 

^o(.«l) = M + L + ^V--^log^- + ^^V- + ... 
4A , oM 4A2 5M A5/2 , gM 

The above result is a string calculation based on the large A assumption. However, as 
discussed in [^ |l^], a general form of E interpolating between weak and strong 
coupling is expected to take the form 

/I T\/f\ 2"+"^ 

S = M + L + L ^ ^c„„(A)A" (-log-) . (2.7) 



n>l m>0 



The coefficients c„_m(-^) have a regular expansion around A = and a strong coupling 
expansion in inverse powers of -s/A. Quite remarkably, some of them are protected and 
are thus genuine constants independent on A. This follows from the comparison between 
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string theory and gauge theory of the 1-loop and 2-loop leading and subleading correc- 
tions [ p8| , |39| , p) , p2[ , |43| ] to the thermodynamical limit of similar circular string solu- 
tions. In particular, this is true for the coefficients 

cioi cii, ci2, C20, andc2i, (2.8) 



and we can write the very explicit expansion (Eq. (1.15) of [p3[]) 



E = M + L[l+ (2.9) 

A , p M , Cii , M C12 , p M , 
^V^(c,o + f log- + ^V- + ...) + 

, .M, C21 , M C22(A) , oM 
+^ log^ ^ (CO + f log - + ^ log^ - + ...) + 

, 6^/ C3l(A) , M C32(A) , 

+^ log^ - (C30(A) + ^ log - + ^ log^ -j; + ■■■)]+■■■, 

A comparison with the one-loop string energy gives (ci_2 would require E2 in the fast long 
string limit) 

1 



cio — ^^20 — --5— J) 

2 TT^ 8 TT* 

(2.10) 

4 4 

3 TT'^ 5 7r° 

The value C3o(0) is discussed in [|^] and is obtained by consistency of the string result 
with the universal dressing phase [p3[]. It should be C3o(0) = 



Now, the crucial point is that we can take the generalized limit Eq. ( |1.4D in the above 
interpolating expansion Eq. ( [2.9[ ). Doing so and assuming the above values for the pro- 
tected coefficients we find the following prediction for the large-j behavior of the FRS 
generalized scaling function 

8 64 1 
j 3 TT p 

512 

7 



/''^(j) = 7-^i + ---. (2-11) 

3 TT 

32 1024 1 

/'^^(;) = ^ + ---. (2.13) 



The two terms in Eq. ( 2.11 ) have actually been already obtained in [gj] by working out the 
finite size corrections to the semiclassical expansion of the sl(2) invariant one-loop spin 
chain. The other two expansions are a higher loop test of the AdS / CFT correspondence 
and have not yet been computed in the gauge theory. We now discuss the confirmation of 
Eqs. ( [2.11| , |2.12D in the context of the large-j FRS integral equation. 

3. The FRS equation in brief 
3,1 Setup 

We consider 5[(2) scaling operators of the form 

O = Tr{D^ Z^) + -- - , (3.1) 
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where D is a specific component of the covariant derivative and Z a scalar field of j\A = 4 
SYM [|^]. The omitted terms are analogous operators with the same number of deriva- 
tives and scalar fields. They are required to form an eigenstate of the dilatation operator. 
As usual, h is also identified with the twist, i.e. the classical dimension minus the Lorentz 
spin, here equal to M. 

The anomalous dimensions of scaling operators of the form Eq. can be orga- 
nized in irreducible multiplets of the st(2) algebra and the top states fill a band, see for 
instance [^^]. We can split the scaling dimension A(5) separating out the classical di- 
mension and define the anomalous dimension 7(5) as 

A(5) = L + M + 7(5). (3.2) 

In terms of the energy E{g) of the s[(2) C psu(2, 2|4) long-range integrable spin chain [Q], 
we have 

7(5) = 2 52£;(5). (3.3) 

The quantity E{^g) is the energy level of an integrable system. Therefore, it is computed 
by solving Bethe Ansatz equations with suitable mode numbers identifying the relevant 
state in the above band. 

In the FRS limit Eq. (^), the Bethe roots u of the minimal state in the band are de- 
scribed by a continuous distribution with density Pm(^i) supported in the region \u\ > 
c(5,j)/ where c is some function of g and j that we shall call gap in the following. The 
label m in Pia{u) stands for magnons, which is the standard name for the excitations of 
the integrable chain. Actually, the relevant quantity in the FRS limit is a specific contribu- 
tion to Pmiu) called (7{u) in and representing a fluctuation component of the magnon 
density. 

Remarkably, there is a dual description in terms of the complementary Bethe roots 
called usually holes. In the FRS limit, the holes are also described by a continuous distri- 
bution with density Ph{u) = p{u) supported in the complementary region |u| < c{g,j). 
The two dual descriptions are fully equivalent and can be connected by the simple relation 

jp{u) = --8a{u). (3.4) 

TT 

The FRS equation is an all-order integral equation for the fluctuation density of magnons 
a{u). It can be turned into an integral equation for the hole density p{u). We shall show 
that this latter equation is better suited for the large-j expansion, at least at the two loop 
level at which we work. 

3.2 The all-loop FRS equation 

We need a few definitions in order to write down the FRS equation. They are fully dis- 
cussed in [^] and we summarize them here for completeness. 
First, we define the BES kernel 

Kit,t') = Koit,t') + Ki{t,t') + Kdit,t'), (3.5) 
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where (Jn(i) is the n-th Bessel function) 

t Ji{t) Mt') - t' Jo{t) Ji{t') 



Ko{t,t') 
K,{t,t') 



t' Ji{t) Jo{t') - t Jo{t) Ji{t') 



l2 



Kd{t, t') = 89^ dt" Ki {t, 2 g t") 



t" 



Koi2gt",t'). 



(3.6) 
(3.7) 
(3.8) 



In this paper, we shall not need the dressing kernel Kd- Then, we define also the hole 
kernel 



Khit,t';c) 



t'-t 
e 2 

t'-t 



du cos{tu) cos{t' u) 



(3.9) 



e 2 t cos{ct') sin(ct) - t' cos{ct) sm{ct') 



2-Kt 



Finally, we define the full kernel /C as 



lC{t,t') = g^K{2gt,2gt') + Kh{t,t';c) 



Jo{2 gt) sin ct' t 



t 2-Kt' 

roo 

-4:g^ / dt"t"K{2gt,2gt")Kh{t",t';c). 



e 2 



(3.10) 



After these preliminary definitions we are ready to write the FRS equation which holds 
for the Fourier transform of the (even) magnon fluctuation density (t{u) 



du cos(< u) a{u). 



The all-loop FRS equation reads 

— 1 



/C(i,0) -4 / dt' JC{t,t') a{t') 



The j parameter is related to the gap parameter c by the relation 

4c 16 /"^ . sin ct 



TT TT Jo 



dt- 



62 a{t). 



(3.11) 



(3.12) 



(3.13) 



The generalized scaling function f{g,j) has a rather complicated all-loop expression 
in terms of the solution to the FRS equation 



f{9,j) = 8g^ 



1-8 l°° dt^^^^^tKh{t,0;c) 



2gt 



2gt 

roo 

a{t) dt'Kh{t,t']c)a{t') 
Jo 



(3.14) 



As shown in pSQ, this can also be written more simply as 

f{g,j)= j + 16 a{0). 



(3.15) 
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4. The two-loop hole density equation in u-space 



The FRS equation can be rewritten in u-space by Fourier analyzing Eq. ( |3.12| ). We did the 
analysis up to the two loop level. After some manipulations we arrive at the following 
result where we use the notation of Appendix for the Ga functions 



2 1 dv 
P(^) = — - :tt:<^i/2(«) + / ir-Go{u-v)p{v) + 

TV J An J — c A 1\ 



(4.1) 



+9' 



TT 



G'i mill) ,^ 

2 7r ^'^^ ' 4; cosh2(7ru) 



71 



In this equation, p(°5 is the one-loop term in the weak coupling expansion of the hole 
density 

p{u) = p^^\u) + g''p^'-\u) + 0{g% (4.2) 
and the fimctional 71 [p] is defined as 

71 [p] = 8 + 2 y ' [Gi/2 {u) + 2 7b] p{u) du. (4.3) 

The generalized scaling function turns out to have the following explicit expression 



fi9,j) = 9^1i[p] + 



(4.4) 



+9' 



8 J -Cs 



TT 



24; 



71 



-Gi'/2(u) + 4C3] p{u)du 



The one-loop terms in Eqs. (|4!l , 4.^ ) are of course identical to those already written in [^]. 
As a check of the two-loops terms, one can compute the small j expansion of Eq. ( ^J| ) and 



reproduces perfectly the results of [ pSQ for the generalized scaling function at two loops. 
Also, the two-loop expression of the gap in this regime is 



c{9,3)=J + 



(4.5) 



+f |^-^log2 + ^(-3 7r3 log2 + 77rC3) | + 
+f log2 2 + ^ (-TT^ + 240 TT^ log2 2 - 672 TT log 2 Cs ) I + 



5. NNLO large-j expansion of the one-loop FRS equation 

We now work out the next-to-next-to-leading (NNLO) large-j expansion of the one-loop 
hole density equation and generalized scaling function. This is feasible and gives the 
value of the leading term in Eq. ( |2.11 ). This is a confirmation of the calculation described 
in [ |2^ ] obtained independently in the large-;' FRS context. Also, we find additional infor- 
mation on the density profile as well as on the dependence of the gap on ; 3> 1- 
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Using the notation of Appendix @, the one-loop hole density p{u) satisfies the equa- 
tion 

p(u) = — --^Gi/2(u) + T^ r dvGo{u-v)p{v), (5.1) 

IT J Z-K ' ZTT J-c 

with the normalization condition relating c and j 

J p{u)du = l. (5.2) 
The one loop contribution to the generalized scaling function is simply (■!/'(z) = ^ log T{z)) 
f^'\j) = 8 + 2; [Gi/2(u) -2tp{l)] p{u). (5.3) 

5.1 Leading order 

If J — > 00 we expect c — > oo. In this limit and using the normalization condition, we can 
write the density equation as 

p{u) = ^ dv [Go{u -v)- Gi/2(u)] p{v). (5.4) 
We now recall the useful integral representation 

^oo^^ e, -cos{ty) cos(ta;) = ;^(Go(a;-y) + Go(a; + y)-2Gi/2(cc)). (5.5) 
Since p{u) = p{—u), we obtain 



oo roo 



p(u) = - I dv I dt ^' cos{tv) (5^) 

TT J-oo Jo e' - 1 

Introducing the Fourier transform 

/oo ^-oo 
due'*"p(u)=/ du cos{tu) p{u), (5.7) 
-oo J — oo 

roo 



1 r°° 1 
p[u) = — dte-'^''p{t) = - I dt cos{tu)p{t), (5.8) 



2 IT J-oo TT Jo 

we find 



This means that 



p(u) = - r dt ^\ ^^^^ cos(i u). (5.9) 
TT Jo e' — 1 



pW = ^V-4^. ^ pW = e-^, (5.10) 
— 1 



where we do not restrict t to be positive. Hence, at leading order. 



The contribution to the scaling function is easily evaluated. We start jErom 



Jo 1 - e ^ 



and compute 



/ dt r-^ — r = / dt T = 0. (5.13) 

This shows that 

/(!)(-,•) = 0- J + (5.14) 

c = 0(j). (5.15) 

and the dominant term linear in j cancels. The constant 8 contribution to f^^\j) is in- 
cluded in the 0(1) terms. 

5.2 Next-to-leading order 

The expansion at large j is not trivial and is better performed in u-space. Let us write 

p{u) = po{u) + 6p{u), po{u) = -3- ^ (5.16) 

The density equation can be written 

2 r dv „ , ^ , ^ r dv 

'j J\v\> 

Rescaling u = cx,v = cy, the first integral can be uniformly expanded as 



2 r dv dv 
5p{u) = — - Go{u - v) po{v) + —Goiu-v)5piv). (5.17) 

7 ^|i;|>c 27r J-c 27r 



where 

#(x) = (1 + ^) log(l + a:)-(l-a:) log(l - a:) ^^^^^ 

2 TT tZ) 

This suggest to set 

Spi-^) = ^ Pi ) + ■ ■ ■ ■ (5-20) 

The density equation is then 

In terms of a;, y it is 

^-^-*(^) + r ^Go[c[x-y))pr[y) = Q. (5.22) 
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Using 



Go{cx) = 2 log c + 2 log |a;| + 
c = a J H , 



with an undetermined constant a, we arrive at the singular problem 

- dy log \x-y\pi{y) = ^{x) , 



dypiiy) 



Notice that the normalization of pi is also consistently predicted by 



(5.23) 
(5.24) 



(5.25) 
(5.26) 



du 



Pq{u) + -^Pi 



+ 



(5.27) 



which, evaluating the elementary integral involving po, reads 



1 = 1 + - 
c 



-- + / dxpi{x) \ + o(^ 

TT j-l J 



(5.28) 



and leads to Eq. (|5.26|) . The constant a can be determined dividing by -s/l — x"^ and inte- 
grating. The key result 



1 log \x - y\ 



-TT log 2, 



together with the elementary integral 



dx 



TT, 



as well as 



leads to 



\x. 



1 log 2 



^/^ 



log 2 1 log 2 



TT 



2 a, 



1 

a = -. 
4 



(5.29) 



(5.30) 



(5.31) 



(5.32) 



TT 2 TT 

The solution of the integral equation with logarithmic kernel is standard. By taking a 
derivative it is reduced to a finite Hilbert transform problem and the solution is 



Evaluating the principal integral we get 



Pi(a;) 



1 + 



1 - 



y-x 



(5.33) 



2 7r Vl - a;2 1 + Vl - a;2' 



(5.34) 
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The evaluation of the generalized scaling function can be done as follows. We want 
to compute 

I du [Gi/2(u) + 27b] (^Po{u) + ^ pi{u/c) + ■ ■ ■ j = (5.35) 

1 

du [Gi/2(w) + 27b] Po{u) + - j dx [Gi/2(cti) + 275] pi(a;) H 



The first integral reads 



<^ ^ ■tp{^ + iu) + ip{l - iu) + 27s 



We know that F(oo) = 0. Also, 

ipil + i + -tc) + 27g _ _2_ 

7r(c2 + |) ~ TTC^^"" " ' ■ - Vc^ 

Hence, 

^■(0) = - — (l+7s + logc) + of4V (5-38) 
The second integral is for large c 

dx [Gi/2(cu) + 27b] Pi (a;) = (5.39) 

2 (7b + log c) J dxpi{x) + 2 J dxlog\x\pi{x) = (5.40) 

2 (7. + log c) ^ f*(0) - J-U ^l2-±Mf) + 1 _ 1. (5.41) 



TT V 2 TT / TT TT 

Combining, we find 

l^^du [Gi/2(ti) + 27b] (^po{u) + ^ pi(«/c) + •••)= + ^ (^) • (^-42) 



We conclude that 

/(^)(j) = 8 + + ^ + =0-j + + oQj , (5.43) 

c = ^ J + subleading. (5.44) 
5.3 Next-to-next-to-leading order 

Expanding further the density equation with the position 

M^^) = ^Pi(f)+^P2(f )+■■-, (5.45) 
and using the results in Appendix, we find 

/ 1 ? 0°^ c + \og\x - y\ + -^6{x - y) + ■ ■ - ^ (^-^ pi{y) + ^ p2(y) + ■ ■ ■ j 
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We assume the following general expansion of the gap 



Hence 



and 



c = ^+;0 1ogj+7 + --- 



j=4c-4;0 1ogc-8;0 1og2-47 + 



+ 



7 /3 log 2 P log c\ 1 



TT j 2lT C \2'K 

The normalization condition 



+ 



+ 



du 



PoM + 7? Pi T + 7TP2 



TT 



U 



27r 



+ 



leads to 



J dxp2{x) = 0. 

The relevant terms in the density equation are thus 

1 f ^ , , I I ^ ^ 1 . ^ 7 P log 2 ;0 log c 
- dx log\x - y\ p2{y) = - pi{x) - — — 

TT J -I 2 2 TT TT 



27r 



(5.47) 

(5.48) 
(5.49) 

(5.50) 
(5.51) 

(5.52) 



The appearance of the logarithmic term log c is tricky and can be understood as follows. 
Dividing the above equation by Vl — a;^ and integrating between — 1 and 1 we find using 
the normalization of p2 



°=2 



1 



^^^^^ ^-Pl0g2-^^. 



dx , 

1 2 



(5.53) 



On the other hand, the pi function has the following behavior for 

1 



Thus the integral in Eq. ( [5.53| ) is singular 



(5.54) 



dx 



Air J -I 1 — x^ 



00. 



(5.55) 



The most singular and universal term is evaluated by integrating over [— 1 + e, 1 — e] and 
identifying e ~ 1/c. The logarithmic singularity is not ambiguous and reads 



^ ^ dx f '^^^^ , = i^i-^ + less singular. 



vT 



x-' 



This gives 



P 



27r 



(5.56) 



(5.57) 
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Apart from this, we shall not attempt to determine more precisely the function p2 nor the 
constant 7 which we shall not need in the end. Instead, we notice the important relation 

i /' <ix log M ^(rt = i ,.(0) - f - ^ - ^ (5.58) 
Computing the central value 

Pi(0) = T^, (5.59) 

4 TT 

we thus obtain a crucial piece in the determination of the anomalous dimension at NNLO 

dx log \y\ p2{y) = ^ - I - ^ log 2 - (5.60) 

If we look back to the expression determining the anomalous dimension, we see that all 
the terms from po are already at the precision of NNLO. Also, the integral involving pi is 

dx [Gi/2 icx) + 2 je] Pi (x) (5.61) 

Due to the results in the Appendix, the expansion of G1/2 has no (5-term and is also already 
at the NNLO precision. The missing piece is (using again the normalization of P2) 

dx [Gi/2(ca;) + 2je] P2{x) = (5.62) 

2 (7b + log c) j dx p2{x) + 2 dx log |a;| p2(a;) = (5.63) 

/"I 1 
2 / da; log I a; I p2 (a;) = 7 - 2/3 log 2 - ;0 log c. (5.64) 



-1 4 
Expanding in powers of j the combination 



we find 



8 + 2j (-- + ^( 4-7- 2^1og 2-/3 logclj+0(^), (5.65) 



f^'\j) = 0-j + + ^ + l^-^y (5.66) 

c = Jj + ^log; + 0(l). (5.67) 
Notice also that Eq. ( 5.66 ) is completely independent on both p and 7 that, to be honest. 



deserve a full determination at next order. 

6. Large-j expansion of the two loops FRS equation: LO and remarks 

The leading order calculation is very simple. We start from the surviving terms for ji — j> 00 
and at this order we can set c = 00. The equation to be solved is 

1 r°° r 1 
p{u) = — dv [Go{u -v)- Gi/2(u) - 5' G'i,^{u)\ p{v). (6.1) 
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A Fourier analysis completely similar to the one discussed for the one-loop case gives the 
interesting result 

The integral of the two-loop correction vanishes and one easily proves the cancellation 
of 0{j) terms in the expression of f^^\j)- The NLO correction is also similar to the one- 
loop case and one proves the cancellation of 0{j) terms in the gap as well as cancellation 
of C(l) terms in f^^\j). We did not push the calculation further since we expect that 
(j) = 0{l/j^) which means that we need a N^LO calculation ! For this reason we shall 
discuss in the next section a fully numerical determination of the two terms in Eq. ( |2.12| ). 
In perspective, this shows that beyond one-loop more sophisticated analytical tools are 
necessary instead of the brute-force one-loop analysis. 

7. Numerical study of the hole density equation 

We look for a numerical determination of the one and two loop densities in Eq. (^). We 
need also to expand the gap 

c = c(0)+52c(i) + O(/). (7.1) 

7.1 One loop 

The numerical problem at one-loop is the solution of the non-singular integral equation 
pW{u) = — -_Gv2(«)+ / ^—Go{u-v)p^°\v)dv, (7.2) 

IT J 2 IT ' J-c(-°^ 2 IT 

1 = / p(°\u)du. (7.3) 

To this aim, we fix c(°) and discretize the u space evaluating the integral by the Boole's 
rule. The integral equation becomes a linear problem which can be solved very efficiently 
and with high accuracy. The resulting density is plugged in the area constraint and c(°) 
is determined by bisection. All the procedure must be repeated with smaller and smaller 
lattice spacings until convergence is achieved. The good convergence is shown in Fig. (||) 
at = 10. 

As a first result, we show in Fig. (Q) the one loop hole density from the numerical 
integration of the FRS equation at the three values of the j parameter j = 10, 30 and 
j = 60. The density is progressively better represented by the LO analytical expression as 
j increases. One notices that the tails of the density, near the boundary of the gap interval, 
show an interesting small rising. This is precisely captured by the NLO solution as shown 
m Fig. d). In that figure, we show the numerical density after subtraction of the analytical 
LO contribution. For the j = 30, 60 we also superimpose the analytical NLO solution. 
It is a rather good approximation, although the finite j data cannot show the divergence 
on the boundary u = ±c. As a further check, we also show, in the j = 30 panel, the 
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numerical solution of the logarithmically singular equation Eq. ( |5.25| ) that determines pi . 
The agreement with Eq. ( |5.34[ ) is of course perfect. 

The dependence of the gap on j is illustrated in Fig. @ where we subtract out the 
leading contribution j'/4 in order to better display the subleading terms. Indeed, a non- 
trivial reminder can be seen which is very well fitted by the heuristic logarithmic term 
Eq. ( |5.57D discussed previously. 

Finally, we show in Fig. the numerical computation of the generalized scaling 
function at one-loop. The data are very well reproduced by the NLO prediction Eq. ( |2.11 ). 
We also show that the LO prediction is not enough to reproduce the numerics. This is 
a confirmation of the NLO contribution. By the way, one can also make a general 2 or 
3-parameter fit to predict a priori the two coefficients and of course they are matched with 
good precision below the 0.1 % level. 

7,2 Two loops 

The two loop equation for p(^) and the constraint on c^-^^ are 

-c(o) 2 7r 2 7r ' 4j cosh (ttu) 

+ J_ p(0)(c(0)) [Go(u - c(°)) + Go(u + c(°))l , (7.4) 

2 TT L -I 

2 p(°) (c(°) ) c(^) + / p(^) (u) = 0. (7.5) 

By discretization, the first equation is a linear problem where we have to insert various 
quantities computed in the solution of the one-loop problem. This must be done with a 
fixed c^^^ which is then evaluated by bisection to impose the second constraint. 

The two-loop contribution to the density profile is illustrated in Fig. (^. Apart from 
remarkably small corrections, the LO expression captures essentially the numerical data. 
The gap is shown in Fig. ^ where we confirm that the two loop contribution starts 0(1/ j). 
We have fitted the numerical results with a 3-parameter fit. The leading term is very 
accurately 6/j. In the figure, we show the result of a 2-parameter fit with the leading term 
fixed at that value. 

Finally, in Fig. (^), we show the two loop generalized scaling function. As in the 
one-loop case, one can predict the coefficients. The result is in perfect agreement with 
Eq. ( |2.12| ). This is best illustrated by superimposing the NLO curve which reproduces 
numerical data very well. Again, one can try to see the accuracy of the LO term alone and 
the figure shows that it is not enough. This means that our computation strongly suggest 
the validity of the expansion Eq. ( |2.12| ). 

8. Conclusions 

In this paper, we have exploited a very simple remark, i.e. the observation that the large- 
j limit of the FRS equation can be used to capture the fast spirming long string limit of 
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AdS^ X superstring perturbation theory. Indeed, in this limit, the string energy can be 
expanded in inverse powers of j as 



E = M + L + /(A, j) log M + C»(M°), (8.1) 
/(A,;) = ^n(A)r". (8.2) 



n>l 



In many respects, the coefficients C„(A) are similar to the more studied scaling function 
and its generalizations which are obtained by expanding f{g,j) around j = 0. Indeed, 
the various C„(A) are defined for all A by the above relation and can be computed at any 
A, weak or strong, by taking the large-j limit of the FRS equation. As we explained, the 
large-j limit is quite natural from the point of view of string perturbation theory where the 
two regimes of slow and fast long string emerges quite symmetrically and are associated 
with j C 1 or j ^ 1. 

As a first attempt to study the large-j limit of the FRS equation, we have described 
in this paper what can be learned in the weakly coupled gauge theory. This regime has 
a non vanishing overlap with the string calculation since part of the C„(A) coefficients is 
protected leading to a prediction from string theory valid also at small A. 

At one-loop in the gauge theory, we can match the result Eq. ( |2.11 ). This is not a new 



check since Eq. ( |2.11 ) has already been obtained by working out the finite size corrections 



to the integrable XXX_i/2 spin chain Nevertheless, this is an important check of 
the approach and several interesting new details are uncovered. In particular, we have 
obtained various results concerning the large-j Bethe roots density and gap dependence. 

At two-loops in the gauge theory, we can match Eq. ( [2.12[ ). This is an interesting check 
first proposed in [^] and never verified. We did not work it out in a fully analytical way, 
but have shown that a numerical approach is feasible and strongly supports a perfect 
agreement. In principle, an analogous study could be carried over to test the three loop 
result Eq. ( ^13] ). 

Clearly, the most interesting development of our analysis is to compute the strong 
coupling expansion of C„(A) from the FRS equation. It would be very interesting to in- 



vestigate whether the effective techniques developed in [ |19t [28[ ] for j <^ 1 can also be 
applied to the large-j case. 
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A. Technical results about the combinations ip{a + ix) + ip{a — ix) 
Let us consider the function itp{z) = ^ log T{z)) 

Ga{x) = tp{a + ix) + tp{a — ix), a, a; G M, and a > 0. (A.l) 
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The function Ga{x) is real. Due to a special reflection property of the ■^(z) function, one 
has the remarkable identity 

Gi{x) = ij{ix) + ij{-ix). (A.2) 



Hence, it will be convenient to extend the definition Eq. ( |A.1| ) to the case a = by under- 
standing 

Goix) = G,{x). (A.3) 

From the integral representation 



coo 



-<P{z) + 1E= ds — —, Re(z) > 0, (A.4) 

Jo 1 - e ^ 

we can obtain the remarkable definite integral 

Gaix) -2 \og\x\ = -2 cosfsa;). (A.5) 

Jo \1 — e ^ s j 

This means that the following Fourier transform holds 

^{G,(a;)-2 1og|a;|} = -27r f -1-!^ - 1 J . (A.6) 



Also, using the above reflection identity, we have for a = 

^{Go(a;)-2 1og|a;|} = -27r (^-J— - 1^ . (A.7) 

If we now want to compute the asymptotic expansion for c — > cxd of the integral 

I = f dxGa{cx) p{x), (A.8) 

we simply add and subtract a logarithm and obtain 

/ = 2 log c J dxp{x) + 2 J dx log \x\ p{x) + J dx {Ga{cx) - log {c\x\)) p{x). (A.9) 

The first term in the asymptotic expansion of the last integral is obtained by writing it as 
the integral of the Fourier transforms of Ga{cx) — log {c\x\) and p and expanding the above 
results. We can compactly write the result as the distributional identity 

Goicx) = 2 1ogc + 2 1og|a;| + -^(a;) + of^y (A.IO) 



c \c 



Gi{cx) = 2 log c + 2 log |a;| + O ( ^ ) . (A.ll) 



The 5-term in the a = case is quite important for the discussion of the main text. 
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Figure 1: Convergence of the one-loop contribution f^^\j) at j = 10 as the number of discretized 
points N in u-space is increased. 
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Figure 2: One loop hole density from the numerical integration of the FRS equation at the two 
values of the j parameter j = 10, 30 and j = 60. In the bottom-right panel, we show a detailed 
view of the case j = 60 where one can appreciate the difference with respect to the lowest order 
density. 
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Figure 3: We show pi which is the one loop hole density from the numerical integration of the FRS 
equation minus the analytical lowest order expression of the density. The remaining curves should 
give at large j the NLO correction computed in the text. For j = 30, we superimpose the analytical 
expression of pi {Hilbert transform label) and the numerical solution of the logarithmically singular 
integral equation that determines it (log. integral equation label). For j = 60 we show the numerical 
data and the Hilbert transform result. 
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10 15 20 25 

j 



Figure 4: One loop gap as a function of j. We show the numerical data minus the leading order 
contribution. We also superimpose a fit with the heuristic logarithmic contribution discussed in 
the text plus an additional subleading constant. 
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Figure 5: One loop generalized scaling function. We show the numerical data and superimpose 
two curves with the NLO and LO analytical expression. 




Figure 6: Two loop hole density from the numerical integration of the FRS equation at the two 
values of the j parameter j = 10 and j = 30. These curves are almost coincident with the leading 
order density. 
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j 



Figure 7: Two loop gap. We show the numerical data and superimpose a reasonable 2-terms fit 
including the leading order whose coefficient is very accurately reproduced by a would-be 3-terms 
fit. 
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Figure 8: Two loop generalized scaling function. We show the numerical data and superimpose 
two curves with the NLO and LO analytical expression. 
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